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We present a systematic study of the effect of Y doping and hydration level on the structural 
transformations of BaCeC>3 based on anelastic spectroscopy experiments. The temperature of the 
intermediate transformation between rhombohedral and orthorhombic Imma phases rises with in- 
creasing the molar fraction lofY roughly as (500 K) x x in the hydrated state, and is depressed of 
more than twice that amount after complete dehydration. This is explained in terms of the effect 
of doping on the average (Ce/Y)-0 and Ba-0 bond lengths, and of lattice relaxation from O vacan- 
cies. The different behavior of the transition to the lower temperature Pnma orthorhombic phase 
is tentatively explained in terms of progressive flattening of the effective shape of the OH - ion and 
ordering of the O vacancies during cooling. 
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I. INTRODUCTION 

The BaCeOs perovskite undergoes three phase trans- 
formations starting from the high temperature cubic (C) 
phase: to rhombohedral (R) at T\ — 1170 K, to or- 
thorhombic Imma (01) at Ti = 670 K and to orthorhom- 
bic Pnma (02) at T3 = 563 K, as determined by neu- 
tron diffraction 1 and by combined differential scanning 
calorimetry, dilatometry and X-ray diffractionP The se- 
quence of phase transformations and the various struc- 
tures are well characterized in the undoped state of 
BaCeOs, and even a quantitative description of the spon- 
taneous strains by means of the Landau expansion of the 
free energy has been presentedP Instead, the situation 
is confused when a trivalent dopant, e.g. Y 3+ , is substi- 
tuted into the Ce 4+ place in order to make the material a 
protonic conductor. The understanding of the influence 
of doping on the phase transitions in the perovskite ionic 
conductors is not only of academic interest. In fact, the 
occurrence of phase transformations, especially if accom- 
panied by ordering of the mobile ionic species, protons 
and O vacancies, is closely related to the mobility of such 
ions and to the durability of the material in applications 
like fuel cells or membranes for gas separation PEI Al- 
though various indications have been reported that the 
transition temperatures i n B aCe03 depend on doping 
and on the hydration state JUSHEI n0 systematic study and 
analysis has appeared yet. According to Raman spec- 
troscopy measurements on variously doped BaCeOs, the 
room temperature structure changes to the more sym- 
metric tetragonal and cubic phases with increasing Nd 3+ 
substitution! 1 -^ but a subsequent neutron diffraction ex- 
periment excludes any significant influence on the room 
temperature orthorhombic structure from Nd dopingP-^ 
For BaCei-zY^Os-s (BCY), a change from the 02 to the 
R structure at room temperature was found at x > 0.2 



with neutron diffraction E21 whereas a later x-ray diffrac- 
tion study did not show such a transition to rhombohe- 
dral at room temperature, but rather impu rity phases 
arising from a more limited Y solubility rangeP^EHl More 
recently, anelastic spectroscopy measurements on BCY 
showed that passing from the hydrated to the outgassed 
state with x — 0.1 lowers the temperature of the Ol-R 
transition by as much as 250 K. 15 Here we present a more 
extensive study of the effect of doping and O vacancies 
(Vo) on the phase transitions in BCY, again based on 
anelastic measurements. An interpretation of the obser- 
vations is proposed, assuming that the driving force for 
the octahedral tilting instabilities is the mismatch be- 
tween too long A-0 and too short B-0 bonds, as usual 
for ABO3 perovskites; a minimal model is adopted for the 
changes with doping and hydration level of the tolerance 
factor and of the lattice relaxation due to V<> 



II. EXPERIMENTAL 

The samples of BaCei_ x Y x 3 with x = 0, 0.02, 0.1, 
0.15, 0.3 were prepared as already described, 16 with start- 
ing powders obtained by auto-combustion synthesis! 1 -^ 
followed by crystallization in air at 1273 K for 5 h. No 
oxide impurity phases were detected by x-ray diffraction 
(XRD) after synthesis for x < 0.15. The sample with 
nominal x — 0.3 was not monophasic, since the solubil- 
ity limit of Y in BCY is lower than 0.3. In fact, impurity 
phases are detected by XRD for x > 0.2 13i and by EXAFS 
for x > 0.17. 14 We did not determine the exact concen- 
tration of Y in solid solution in the x — 0.3 sample and in 
what follows we will set this value to x = 0.2. The pow- 
ders were first uniaxially pressed at 50 MPa and then 
isostatically pressed at 200 MPa obtaining 60 x 7 x 6 mm 
ingots, which were sintered at 1773 K for 10 h. The 



samples were cut as thin bars about 1 mm thick and 
~ 4 cm long. In order to make them conducting for the 
anelastic experiments, their faces where covered with Ag 
paint, or Pt paint when temperatures higher than 900 K 
had to be reached. We tried with SPI-Chem Conductive 
Platinum Paint consolidated at 1270 K, or with 1000 A 
of Pt magnetron sputtered directly on the sample sur- 
face. Unfortunately, none of these electrodes resisted the 
anelastic measurements in vacuum < 10~ 5 mbar in the 
temperature range 1000-1300 K, since in all cases they 
evaporated away. This fact rendered the measurements 
at > 1000 K difficult, and we could not obtain extensive 
and reliable data on the R-C transition near 1200 K. 

Hydration was achieved by maintaining the samples for 
1-2 h at 793 K in a static atmosphere of 50 — 100 mbar 
H 2 0, followed by slow cooling, while outgassing was 
achieved in vacuum < 10 -5 mbar up to 1000 K or dur- 
ing the anelastic experiments. The resulting variations 
of the gaseous contents were monitored from the change 
of weight. The reaction of equilibrium of the O defi- 
cient perovskite with water vapor is H^O+Vq'+Oq < — > 
20Hq, 4 where a molecule of water fills one Vo and 
provides two protons that may diffuse among 2 ~, to 
which are bound as peroxide ions (OH) - . The super- 
script dots represent excess +e charges of the species 
with respect to the perfect lattice in Kroger- Vink nota- 
tion. According to this reaction, the concentration of Vo 
in BaCei-zY^Oa-^Hy can vary within 6 < x/2 and con- 
comitantly the concentration y of protons within y < x. 
It was found that the maximum possible hydration was 
~ 15% lower than the theoretical maximum of y — x; this 
is usual in the doped cerate and zirconate perovskites, 
and can be due to partial occupation of the Ba sites by 
the trivalent dopants or to other defects. 

The Young's modulus E was measured by electrostat- 
ically exciting the flexural modes of the bars suspended 
in vacuum on thin thermocouple wires in correspondence 
with the nodal lines. Besides the 1st flexural mode, the 
3rd and sometimes the 5th modes, with frequencies 5.4 
and 13.3 times higher, could be measured during the 
same run; the frequency of the fundamental mode was 
uj/2-k ~ 1.5 — 3 kHz, depending on the sample shape and 
state. The data will be presented as real part s' of the 
elastic compliance s (w,T) — s' — is" = E^ 1 , referred to 
its extrapolation Sq to K, and elastic energy loss coeffi- 
cient Q^ 1 = s" / s' . The first is proportional to the square 
of the sample resonance frequency, s' (T) ex w 2 (T) , and 
presents peaks or steps at the structural phase trans- 
formations; the latter was measured from the decay of 
the free oscillations or from the width of the resonance 
peak, and presents peaks due to the relaxational motion 
of point and extended defect*^ (Vo, protons and their 
complexes with dopants, twin walls, etc.). For an ele- 
mentary relaxation process it is^ 
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with maximum at the temperatures at which the defect 
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FIG. 1: Real part of the elastic compliance (upper panel) 
and elastic energy loss coefficient (lower panel) of undoped 
BaCeC>3, measured at 1.6 and 8.4 kHz. 



relaxation time r ~ ut^ 1 . Since r (T) is a decreasing func- 
tion of temperature, usually according to the Arrhcnius 
law, the temperature of a thermally activated peak in- 
creases with frequency. The peaks due to the hopping of 
Vo and to the reorientation of pro tons around Y dopants 
have already been identified) 15 * 16 ^ and allow one to moni- 
tor the concentrations of such defects and to study their 
dynamics. 



III. RESULTS 

Figure [T] presents the anelastic spectrum of a sample 
of undoped BaCeOs measured at two frequencies: 1.6 
and 8.4 kHz. The real parts s' are practically coinci- 
dent at both frequencies and present sharp steps or peaks 
in correspondence with the three phase transformations 
at Ti = 1210 K, T 2 = 668 K and T 3 = 540 K. These 
temperatures are close to those determined by neutron 
diffraction^ and define the temperature ranges of the cu- 
bic, rhombohedral and two orthorhombic phases. We 
identify T± of the R-C transformation with the temper- 
ature of the kink between almost flat and sharply rising 
compliance, rather than with the peak at a temperature 
65 K lower. This also coincides with a sharp dip in the 
absorption. Such a dip is rather anomalous, since usu- 
ally one finds a peak or more or less rounded step at the 
onset of a structural transformation, but it clearly sep- 
arates the rhombohedral region, with absorption due to 
the motion of domain walls, from the cubic region with- 
out appreciable anelastic losses. 

The elastic energy loss coefficient, besides clear anoma- 
lies in correspondence with the transitions, has five relax- 
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FIG. 2: Real part of the elastic compliance (upper 
panel) and elastic energy loss coefficient (lower panel) of 
BaCei-^Y^Oa-^Hy with x — 0.1, measured at ~ 1.5 kHz at 
various hydration levels, from fully hydrated (black) to fully 
outgassed (light grey). 



ation peaks in the 300 — 1300 K temperature range, la- 
beled R4-R8 because there are other relaxation processes 
at lower temperatures (see Ref. [TBI and Fig. El later on). 
The thermally activated character of these processes is 
clear from the fact that they are shifted to higher tem- 
perature at the higher frequency (see Eq. (fll)). 

Figure [2] presents a series of anelastic spectra of a sam- 
ple of BaCei-zYzOa with x = 0.10 at various stages of 
hydration, from fully hydrated (thick black lines) to fully 
outgassed (light grey); these spectra have already been 
published in a preliminary studjff^ of the effect of vary- 
ing hydration on the structural and elastic properties of 
BCY. The two Q _1 (T) peaks at lower temperature are 
labeled as RH and RV, since they are due to hopping 
of protons, likely around Y dopants™ and of Vo, re- 
spectively. Their evolution allows us to confirm that the 
sample passes from fully hydrated (RV is absent and RH 
saturated) to fully outgassed (RH is absent and RV satu- 
rated). The presence of the intermediate curves (only few 
of them are reported) allows us to ascertain that indeed 
outgassing shifts the transition at T 2 of 250 K to lower 
temperature, while the transition at T3 is soon smeared 
and masked by the presence of the former transition, but 
does not seem to shift appreciably. 

The effect of Y 3+ doping on the phase transformations 
at T2 and T3 is shown in Fig. [3J where the elastic com- 
pliance curves are plotted of samples having x = 0, 0.02, 
0.10, 0.15 and 0.2 in the fully hydrated and fully out- 
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FIG. 3: Real part of the elastic compliance of a series of 
samples with different doping x measured in the hydrated 
and outgassed states. 



gassed states. There is no outgassed curve at X = 0.2, 
because the sample broke after the first measurement. 
With increasing doping, and hence lattice disorder, there 
is progressive smearing of the peaks at the transitions, so 
that T 3 becomes more and more difficult to determine. 
This is especially true in the outgassed state, where the 
transition at T2, whose effects on the elastic compliance 
are prevalent, shifts consistently to lower temperature 
and masks the effects of the 02-01 transition. Note that 
there is an inversion in the trend of the anelastic spec- 
tra between x = 0.1 and 0.15, as discussed later. The 
transition temperatures deduced from these curves will 
be plotted in Fig. [5j 

We finally present an example of <5 _1 (T) curves mea- 
sured at three different frequencies, where it is partic- 
ularly clear that the dynamics of the Vo is not simply 
that of independent defects, which would give rise to a 
Debye relaxation, Eq. (fTl), but seems to have an impor- 
tant contribution from cooperative effects, possibly con- 
nected with Vo ordering in the 02 phase. The Q^ 1 (w, T) 
curves in Fig. [4] are measured at 3.3, 18 and 44 kHz on 
a sample with x = 0.02 in an intermediate state of hy- 
dration where both Vo and H are present. Among the 
various peaks, with the help of Fig. [2] we recognize RH, 
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FIG. 4: Elastic energy loss coefficient of BCY with x = 0.02 in 
an intermediate hydration state, measured at three frequen- 



probably due to reorientation of H about Y dopantsp^ 
while RV is certainly due to Vo, although it is not yet 
determined whether trapping by Y has a role. The in- 
teresting feature in Fig. [4] is that the intensity of peak 
RV is a drastically increasing function of temperature, 
instead of having the \/T dependence expected from Eq. 
(fTl). The effect is not due to O loss during the measure- 
ments in vacuum, since the curves at different frequencies 
are measured during a same run, and we have abundant 
data showing that the Q^ 1 (T) are perfectly reproducible 
until one does not exceed 500 K. We mention that also 
the relaxation peak R6 in Fig. [2] displays a similar be- 
havior, although the divergence of its intensity on ap- 
proaching T3 might be partially due to overlapping with 
the narrow dissipation peak associated with the struc- 
tural transition. We did not make a thorough analysis of 
these complicated anelastic spectra, and cannot say yet 
whether R6 involves Vo, twin walls or both; it certainly 
appears that the motion of the Vo has a high degree of 
cooperativity below T 3 . 



IV. DISCUSSION 

Figure [5] shows the transition temperatures T\ , T 2 and 
T 3 plotted versus Y doping in both the fully hydrated and 
outgassed states. The temperatures T 2 and T 3 are deter- 
mined from both the real parts s' (T) in Fig. p] and the 
respective Q~ x (T) curves (not reported here), while T\ is 
determined from Fig. [l] The Ol-R phase transformation 
is the better characterized in the present measurements 
and exhibits the largest dependence on doping: an in- 
crease of T 2 with doping in the hydrated state and an 
even larger decrease in the outgassed state. Notice that 
there is a difference of nearly 30 K between outgassed 
and hydrated state even at x = 0, which may be due 
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FIG. 5: Temperatures of the three structural transformations 
of BCY in the hydrated (filled symbols) and outgassed (open 
symbols) states: R-C square, Ol-R triangles, 02-01 circles. 
The thick filled and dashed lines are calculated as explained 
in the text. 



to the contribution of electronic compensation or to the 
presence of some defects, e.g. from non perfect stoichiom- 
etry. The transition at T3 between the two orthorhombic 
structures is less affected by doping and by O stoichiom- 
etry. Differently from that at T 2 , it decreases slightly 
its temperature with doping in the hydrated state, while 
only at x = 0.02 it was possible to verify the lowering of 
the transition temperature after outgassing, because of 
the masking effect of the transition at T 2 . 

The R-C transformation has only one point for the un- 
doped case, because our results are only partial and pre- 
liminary, due to the experimental difficulties explained 
in Section II. In addition, our anelastic experiments are 
made in high vacuum, so that above 800 K it is impossi- 
ble to maintain the sample in the hydrated state, and we 
can only measure reliably the temperature T\ of BCY in 
the outgassed state. Also in the literature there are no 
data on T x of doped BCY. 

In what follows we will try to explain the fact that 
the temperature T 2 (x) increases with Y doping x in the 
hydrated state, and instead decreases with x of an even 
larger amount in the outgassed state. The uncertainty in 
the values of T 2 in Fig. [5] is smaller than the symbol size, 
and the fact that the points at x = 0.10 and 0.15 do not 
follow a monotonic trend with doping is likely real and 



not an experimental vagary. This anomaly can be put in 
relation with the observation of a jump or extremum in 
the doping dependence of several structural parameters 
of BaCe03_a and SrCeC^-^ at a nominal concentration 
of Vo <5 ~ 1/16,121 corresponding to 6 = 2x = 1/8 in the 
fully outgassed state. This phenomenon has been ten- 
tatively attributed to ordering of the Vo commensurate 
with the lattice during the synthesis at high tempera- 
ture, hence with possible ordering of the cation dopants 
that would affect the hydration properties also at lower 
temperatureP21 We will ignore this local inversion of the 
variation of Ti [x ~ 1/8) , and only consider the positive 
average derivative of T^ (x, S — 0) with respect to x, and 
negative derivative of Ti (x, 8 = x/2). 

In searching for the relevant factors determining Ti , we 
note that the Ol-R transition involves tilting of the O oc- 
tahedra, without the atomic off-centering accompanying 
the ferroelectric transitions or additional Jahn- Teller dis- 
tortions, since neither Ce 4+ nor Y 3+ are Jahn- Teller ac- 
tive. BCY is also inert from the magnetic point of view, 
so that we conclude that the main driving force for the 
octahedra to tilt is the mismatch between too long B-0 
bonds (B = Ce/Y) and too short A-0 bonds (A = Ba), 
as usual for AB0 3 perovskitesPS The tendency of per- 
ovskites to undergo tilting transitions is often expressed 
in terms of the tolerance factor 



TABLE I: Ionic species of BCY and their molar fractions, 
radii and coordination numbers, according to Shannon 
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where the mean ionic radii are the effective ones usually 
taken from Shannon's tables P3 A value t — 1 means that 
the ideal A-0 and B-0 bond lengths, taken as the sums 
of the ideal ionic radii, perfectly match the cubic struc- 
ture, and therefore that the cubic phase should be stable; 
t < 1 means that the B-0 bond length is too large with 
respect to the A-0 one, and therefore that the octahedra 
tend to rotate in order to accommodate the mismatch. 
The longer and weaker A-0 bonds have larger thermal 
expansion than the shorter and stronger B-0 bonds. For 
this reason, perovskites with t < 1 already at high tem- 
perature further decrease t on cooling, until the mismatch 
between too long B-0 and too short A-0 bonds is relived 
by a tilting structural transformation. Usually, with de- 
creasing t below 1, one finds first tilt patterns producing 
a more symmetric rhombohedral struct ure an d then the 
more distorted orthorhombic structures ! 19 * 21 ! In this re- 
spect, BaCe03 behaves normally, with the cubic phase 
transforming into rhombohedral R3c and further into or- 
thorhombic 01 (Imma) and 02 (Pnma). The tilt pat- 
terns in Glazer's notatiorpJ^SJ are respectively a~a~a~, 
a°b~b~ and a + b~b~ , and the anomaly in the sequence of 
transformations is the intermediate loss of a tilt system 
passing from a~a~a~ to a?b~b~ . Yet, the general trend 
of C, R and O structures with decreasing t is obeyed, and 
the final a + b~b~ tilt system is the usual ground state of 
tilted perovskitespS also favor ed by the slightly covalent 
component of the A-0 bondsP^EI] it can be concluded 
that the tolerance factor should be the relevant parame- 
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tcr in promoting the structural transformations in BCY. 
Another indication in this sense is the fact that SrCe03, 
having a still smaller t due the smaller Sr ionic radius, 
remains in the 02 phase at least up to 1270 K. 26 

It has been noted that, in perovskites with cation 
chemical disorder in the A sublattice, the temperatures 
of the structural, and especially magnetic and electronic 
transitions appear to be sensitive to both the tolerance 
factor, which measures the coherent strain effect, and the 
variance of the A cations sizes, which measures the in- 
coherent partpB In the present case we are dealing only 
with structural transformations, without the additional 
critical dependence on the bond angles involved in the 
electronic and magnetic transitions, and we will just take 
into account the average effects included in t. 

For the limited objective of understanding the effect of 
doping and hydration on the structural transformations, 
but not their detailed nature, we keep the analysis as sim- 
ple as possible, following the idea that the driving force 
for the n-th tilting transition is the decrease of t below 
some critical value t n , and therefore that the transition 
temperatures are proportional to such a driving force, 



T n (x) = ATx(t n -t (x)) 



(3) 



The dependence of t on x can be estimated by assum- 
ing Vegard's law, namely that the introduction of a mo- 
lar concentration x of defects, each contributing with a 
change Sv to the ionic volume, causes an isotropic vol- 
ume change equal to AV — 5v x. In addition we consider 
a purely ionic picture with each ion having its nominal 
valence. The ionic species and their molar fractions, co- 
ordination numbers and radii are listed in Table I. 

Notice that there is no difference between the use of 
ionic and crystal radii of the Shannon t ablest since they 
all differ by ±0.14 A, depending whether they are an- 
ions or cations, and the A-0 and B-0 ideal distances 
are unaffected by the choice. The tolerance factor of 
undoped BaCeOa resulting from Table I, to be consid- 
ered as referred to the 02 room temperature structure, is 
to = 0.9428. Dealing with complete outgassing or hydra- 
tion, we discard the electronic compensation and assume 
that the chemical formula of BCY with 5 Vo and 2-H2O 
is Ba 2+ Ce 4 + a .Y 3+ 2 I i5 _ OHy , with the charge compen- 



sation requiring 
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Hydrated and outgassed states are therefore defined by 
y = x, 5 = and y = 0, 5 = a;/2, respectively. The pro- 
ton is assumed to form the hydroxide complex (OH) - , 
whose radius is also tabulated. The assumption is cor- 
roborated by the observation that H fills the hole at the 
top of the bonding Ce 4 /— O 2p valence band, mostly of 
O 2p character, which is introduced by trivalent doping 
in dry atmospherePZl Certainly, the approximation of a 
spherical (OH) - ion is inadequate, but this is discussed 
later on. 

The tolerance factor of doped BCY can then be written 
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where Are = ry — rc > 0, Aro = »'oh — r o < is 
not very influent because it appears in the same manner 
in the numerator and denominator of i, and the pres- 
ence of Vo is not taken into account yet. It is some- 
times assumed that the radius of an Vo in a perovskite 
is the same as that of the O 2- ionpSl and the fact that 
O deficient perovskites such as LaCo03_^ increase their 
volume with increasing 6 is attributed to the enhanced 
radius of the reduced B cationPSJEED j n the present case, 
it is evident that a similar assumption would not explain 
the marked depression of the transition temperatures of 
the outgassed state with respect to the hydrated and the 
undoped states. Such a depression would require an in- 
crease of t in Eq. (pi), which is not supported by any 
indication. The reduction of T n must arise from the elim- 
ination of the B-O-B and A-O-A bonds, whose rigid net- 
works compete against each other, with expansive and 
compressive pressures respectively. The introduction of 
Vo therefore relaxes the driving force for the tilting struc- 
tural transformations, reducing it by an amount R (<5). In 
the absence of a more detailed and quantitative estimate 
of the structural relaxation introduced by Vo, we will 
assume that R (S) = (1 — / <5), resulting in 
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where / is a parameter that quantifies the amount of 
lattice relaxation associated with Vo; / = 1/3 would 
correspond to a situation in which a tilting driving force 
exists even with few sparse bonds, and therefore it must 
be / ;g> 1/3. The value of /, or more properly the shape 
of the function R (5) , and particularly the value of 6 at 
which it vanishes, should be connected with the critical 



concentration of Vo at which continuum sequences of 
bonds disappear over some length scale. Yet, there are 
many factors involved, for example at S — 0.5 there might 
be ordering into the brownmillerite structure) 31 * 32 * and 
vanishing of R (5) at 5 = 0.5 would correspond to / = 2. 
We will leave / as a free parameter whose value must be 
> 1/3. 

From Eq. Q we may analyze the various factors pro- 
ducing a variation of T n with doping and hydration. In 
the hydrated state R = 1 and the relevant quantity is, to 
first order in the changes of the ionic radii, 
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where the first positive term, representing the aver- 
age increase of the B radius on doping, is dominant. 
The second term is reduced by a geometrical factor 
~ | (l — l/\/2) ~ 0.1, and therefore the most question- 
able assumption of adopting the tabulated radius for the 
hydroxide ion is not important. From Eqs. (El) and ([5]) we 
obtain the proportionality factor between tilting driving 
force and structural transition temperature, as 



AT 



dT n I 
dx I hydr 



d(t n —t) 
dx 



hydr 



Setting 



^H , = 520 K we obtain AT = 44300 K and 

dx I hydr 

extract the critical tolerance factor t% for the Ol-R tran- 
sition at T2 = 660 K in the undoped case from Eq. Q 
as ti = 0.9577. Finally, the parameter / is deduced from 
the initial slope of ^1 .... '- -1100 K as 
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The resulting T2 (x) curves in the hydrated and outgassed 
states are plotted as thick solid lines in Fig. [5] There 
is some arbitrariness in the choice of the initial slopes 
of the curves, due to the above mentioned anomaly be- 
tween x = 0.1 and 0.15, but the main features can be 
reproduced with reasonable parameters. This simple rea- 
soning might be applied also to the other two transitions 
at T\ and I3, with different values for the critical toler- 
ance factors t\ > t<z > £3 and possibly also for the pro- 
portionality factor AT, because the different structures 
relax the mismatch between A-0 and B-0 sublattices at 
varying degrees. The dashed lines in Fig. [5] are obtained 
letting AT and / unchanged and setting t\ = 0.9701 in 
order to reproduce T\ (0) = 1210 K. It is reassuring to 
find that t\ is exactly at the lower limit of the range 
0.97 < t < 1 where cubic perovskites are found.- 3 - 3 - There 
are no data for T\ in the hydrated state, while for the out- 
gassed state there are only preliminary anelastic spectra 
suggesting that at x = 0.15 it is T x > 950 K, about 100 K 
higher than the dashed line. We refrain from speculating 
whether this would be due to a larger value of AT for 



that transition or to other reasons that are not included 
in the present minimal model. 

Additional factors are likely present in the 01-02 tran- 
sition at T3, whose temperature even decreases slightly 
on doping, maintaining the anomaly between x = 0.10 
and 0.15. Yet, the effect of Vo is again to depress the 
transition temperature, although this is verifiable only at 
x = 0.02, due to the overlapping with the Ol-R trans- 
formation. Among the phenomena interfering with the 
01-02 transition is the ordering of the Vo- In fact, while 
the Vo are disordered in their three equivalent O sublat- 
tices in the rhombohedral structure and also in the or- 
thorhombic 01, they are confined to only two sublattices, 
avoiding the third crystallographically inequivalent sub- 
lattice, in the 02 structurePJ^It is not clear whether Vo 
ordering is concomitant with the transition or it occurs 
at a slower rate after the transition is completed. We 
have already noteoP^ that a possible sign of cooperative 
ordering of Vo is the anelastic relaxation process labeled 
as R6, whose intensity seems to diverge on approach- 
ing T 2 from below, as expected from critical ordering of 
the elastic quadrupoles associated with the VoP^"^ A 
similar divergence of the relaxation strength A is shown 
in Fig. [4] for peak RV. In the framework of the Bragg- 
Williams approximation} 3 ^ the critical temperatur e for 
the onset of ordering of Vo is expected to scale a d 35 * 36 ! 
6(1 — 6) and therefore to increase with doping, possi- 
bly driving the 02-01 transformation to higher temper- 
ature and explaining the apparently reduced decrease of 
T 3 (x, S — x/2) with respect to T 3 (x, S — 0), compared to 
T 2 . 

The fact that the decrease of T3 in the hydrated state is 
smaller then for T 2 , instead, must involve completely dif- 
ferent mechanisms. A possibility is that on cooling the 
proton localizes itself more and m ore wi thin the plane 
perpendicular to the B-O-B b on( jUJ12138| effectively re- 
sulting in an increased flattening of the hydroxide ion 
and hence in a reduction of its effective radius along the 
B-O-B bond direction. Such an effect, namely an addi- 
tional reduction of the B-OH-B but not of the A-OH-A 
bond lengths on cooling, would reduce the mismatch be- 
tween the two bond networks and result in a stabilization 
of the higher temperature phase, hence a decrease of T 3 . 

Finally, let us compare these dependencies of the tran- 
sition temperatures on Y doping with those measured by 
XRD and dilatometry with Yb doping. 9 Those data have 
been considered insufficient, 2 due to the smallness of the 
anomalies in the linear expansion, which also exhibit an 
additional dip not associated with any phase transfor- 
mation, and the limited number of diffraction peaks that 
were analyzed. Yet, Yamaguchi and YamadaP plotted 
Ti, T 2 and T3 versus Yb doping measured both under 
wet and dry conditions, as in our Fig. [5] Similarly to the 
present results, the T„ in dry atmosphere are lower than 
those under in wet atmosphere, especially for x > 0.1 
but to a lesser extent than in Fig. [5] The main differ- 
ence between the two sets of experiments is that T 2 (x) 
with Yb has a much weaker rise with doping than with 



Y and only for x > 0.1, while T\ (x) even decreases with 
doping. This difference can also be explained within the 
above model, since the radius of Yb 3+ is slightly smaller 
than that of Ce 4+ , instead of larger as for Y 3+ , so that 
with Yb it is Ar B = -0.003 A instead of +0.03 A. 



CONCLUSIONS 



The temperatures T„ of the structural transforma- 
tions in BCY have been systematically deduced from the 
anelastic spectra as a function of Y doping x and hydra- 
tion level y, or O deficiency 5. The most complete data 
are for the intermediate transition between orthorhombic 
Imma and rhombohedral at T 2 ; the data of the transition 
at the lowest temperature T3 to the orthorhombic Pnma 
phase are incomplete in the outgassed state, due to over- 
lapping with the transition at T 2 , whose effects prevail in 
the anelastic spectra. Of the transition to the cubic phase 
we could measure only the temperature T\ = 1210 K in 
the undoped state. 

The main result is that T 2 (x) increases roughly as 
(500 K) x x in the fully hydrated state and decreases twice 
as much in the fully outgassed state. An anomaly with re- 
spect to the monotonic trend between x = 0.10 and 0.15 
is associated with similar anomalies already observed in 
various structural parameters at the same doping level, 
while the average trend is explained with a simple model. 
As usual, it is assumed that the main driving force for 
the structural transformations, which consist of rotations 
of the B06 octahedra (B = Ce,Y), is the mismatch be- 
tween the more rigid and compressed network of B-O-B 
bonds and the network of Ba-O-Ba bonds under expan- 
sion. The transition temperatures T n are assumed to be 
proportional to (t n — t) (1 — / 8), where t is the tolerance 
factor measuring the ratio between the ideal Ba-0 and 
B-0 bond lengths, t n a critical value of t below which 
the n-th tilting transition occurs, and / a parameter de- 
termined by how much the mismatch stress between dif- 
ferent types of bonds is relieved by the presence of O 
vacancies. In this manner it is possible to explain the 
experimental T 2 (x, 6) data with reasonable parameters, 
and also to reproduce T\ (x = 0) with i x = 0.97, which is 
just the lower limit of the known range 0.97 < t < 1 for 
cubic perovskites. 

The transition temperature T3 , instead, decreases with 
doping and has a reduced difference between hydrated 
and outgassed states. These differences with respect to 
T 2 (x, 6) are tentatively explained in terms of a reduc- 
tion of the mismatch between the bond lengths, due to 
a flattening of the effective shape of the hydroxide ion 
perpendicularly to the B-O-B bond during cooling, and 
to the ordering of the O vacancies in the Pnma phase. 
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